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ABSTRACT 
The  statistical  theory  of  turbulence  is  formulated  in  terms  of  the 
characteristic  functional  r  of  a  probability  distribution  of  velocity  fields 
in  space-time.   A  functional  differential  equation  for  r  is  derived  from  the 
Navier-Stokes  equation  and  the  condition  of  incompressibility.   All  the  moments 
of  the  velocity  field  may  be  determined  as  functional  derivatives  of  r.   The 
n-point  space-time  distribution  functions  of  the  velocity  field  (n=l,2, ... ) 
are  directly  determined  from  the  values  of  T   for  particular  argument- functions. 
The  present  theory  generalizes  that  of  E.  Hopf  in  that  non- simultaneous  as  well 
as  simultaneous  correlation  functions  are  obtained  from  I\   In  addition, 
external  driving  forces  are  included  in  the  formalism.   This  permits  the 
representation,  as  functional  derivatives  of  r,  of  the  Green's  functions  which 
describe  the  averaged  response  of  the  velocity  field  to  infinitesimal  perturbatior 
As  an  illustration  of  the  formalism,  an  explicit  solution  of  the  functional 
equation  for  r  is  obtained  in  the  limiting  case  of  vanishing  Reynolds  number, 
and  the  known  results  of  the  theory  of  "weak  turbulence"  are  thereby  recovered. 


1.  Introduction  and  Summary 

A  functional  formalism  for  turbulence  which  evolves  in  time  according 

r2  3i 

to  the  Navier-Stokes  equation  has  been  given  by  Hopf,  L  '   J  who  introduced 
the  characteristic  functional  $  of  a  probability  distribution  of  simultaneous 
velocity  amplitudes  and  derived  a  functional  differential  equation  for  $. 
If  the  solution  of  this  equation  can  be  determined,  one  can  obtain  from  it 
all  correlation  functions (moments  )  of  the  velocity  field  which  involve 
simultaneous  time -arguments. 

In  the  present  paper,  we  give  a  generalization  of  Hopf's  formalism 
in  which  the  central  role  is  played  by  a  characteristic  functional  T   of  the 
full  space-time  distribution  of  the  velocity  amplitudes.   This  permits  the 
determination  of  the  non -simultaneous  as  well  as  the  simultaneous  correlation 
functions.   As  a  further  generalization,  we  include  driving  forces  in  the 
Navier-Stokes  equation.   In  Section  2,  we  define  T   and  show  how  the  correlation 
functions  are  obtained  from  r  by  functional  differentiation.   We  also  show 
how  the  n-point  space -time  probability  distributions  of  the  velocity  field 
(n=l,2,... )  can  be  expressed  in  terms  of  T. 

In  Section  3>  we  derive  a  functional  differential  equation  for  V   from 
the  Navier-Stokes  equation  and  the  incompressibility  condition.   The  formalism 
is  presented  first  in  terms  of  a  space -time  representation  and  then  is 
transformed,  in  Section  h,    into  a  wave  vector -time  representation.   As  an 
unexpected  by-product  of  our  attempt  to  generalize  Hopf's  work,  we  find  that 
the  new  formalism  is  somewhat  simpler,  in  that  certain  complications  associated 
with  the  elimination  of  the  pressure  term  from  the  equations  are  avoided. 

The  application  of  the  formalism  is  illustrated  in  Section  5,  where 


we  obtain  an  explicit  solution  of  the  functional  differential  equation  in 
the  limiting  case  where  the  characteristic  Reynolds  number  of  the  turbulence 
vanishes  and  the  Navier-Stokes  equation  becomes  linear.   The  solution  is 
obtained  by  analogy  to  elementary  methods  of  solving  linear  ordinary 
differential  equations.   From  our  solution  we  recover  the  known  results  of 
the  theory  of  "weak  turbulence". 

In  Section  6,  we  express,  in  terms  of  r,  a  sequence  of  Green's 
functions  which  describe  statistically  the  response  of  the  velocity  field 
to  infinitesimal  perturbations  in  the  driving  forces.   As  an  illustration, 
we  evaluate  the  Green's  functions  in  the  zero  Reynolds  number  limit. 

The  analysis  outlined  above  is  restricted  to  the  case  of  non-random 
driving  forces.   The  extension  to  stochastic  driving  forces  is  given  in 
Section  f.      We  conclude  the  paper  with  a  brief  description  of  how  the 
conditions  of  statistical  homogeneity  appear  in  our  formalism  (Section  8). 

The  motivation  of  the  present  paper  is  essentially  that  of  Hopf's 
original  paper:  It  is  hoped  that  the  compact  expression  for  the  dynamical 
equations  which  the  functional  formalism  affords  may  suggest  novel  approximation 
schemes  that  would  not  be  apparent  from  the  equivalent  infinite  set  of 
moment  equations.   So  far  this  hope  has  not  been  realized.   The  obvious 
procedures  suggested  by  the  functional  formalism  are  expansion  in  moments 
and  expansion  in  cumulants.   Both  of  these  procedures  are  also  obvious  from 
a  moment  formulation.   We  hope,  however,  that  the  present  enlargement  of  the 
functional  formalism  to  include  the  full  space-time  distribution  and  the 
Green's  functions  of  the  velocity  field  may  provide  a  context  in  which  new 
approximations  can  be  developed. 
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2.  Characteristic  functional  and  space-time  correlation  functions 

The  statistical  theory  of  turbulence  may  be  described  in  terms  of 
a  probability  distribution  P  in  the  space  of  velocity  vector  fields 


(l)       u  =  (u  ,u  ,u  )   =  u(t,x  ,x  ,x  )     =  u(t,x) 


which  satisfy  the  system  of  equations 


(2)  vf  +  Pa     =     -uV     +vua         +fa;  a-1,2,3   ; 

X  XXX 


(3)  uaa      =     0. 

X 

Here  p  is  the  scalar  pressure,  v   is  the  kinematic  viscosity  coefficient, 

GO,     N 

f  (t,x)  is  an  arbitrary  prescribed  driving  force  which  is  defined  for  t  >  0, 
and  the  density  is  assumed  to  be  unity.   We  employ  the  summation  convention 
over  repeated  indices.   (2)  is  the  system  of  (forced)  Navier -Stokes  equations, 
and  (3)  is  the  condition  of  incompressibility.   We  assume  that  the  solution 
u  is  uniquely  determined  for  t  >  0  by  specifying  initial  conditions  at  t  =  0, 
and  appropriate  boundary  conditions  at  spatial  infinity.   Unless  explicitly 

indicated,  integrations  with  respect  to  t  will  be  over  the  semi -infinite 

12?  12? 

interval  0  <  t.   Integrations  with  respect  to  x  =  (x  ,x  ,x  )  and  k  =  (k  ,k  ,k  ) 

will  always  be  over  the  entire  3 -dimensional  space. 

For  functionals  F[u]  (which  assign  a  numerical  value  to  each  vector 
field  u)  we  define  the  mean  value 


CO       <F>  =   [  F  [u]  dP  . 


-  k 


The  integration  here  is  over  the  entire  space  of  velocity  vector 

fields  which  satisfy  (2),  (3)  and  the  appropriate  boundary  condition.   The 

probability  is  non-negative  and  normalized,  i.e.,   /  dP  is  non-negative  if 

A 
A  is  any  subset  of  the  space,  and  is  equal  to  one  if  A  is  the  whole  space. 

It  follows  that  the  probability  distribution  cannot  assign  a  negative 

probability  to  any  velocity  field. 

Now  let 


(5)  y  =  (y\y2,y3)  =  y(t,x; 


be  an  arbitrary  real  vector  field  which  vanishes  at  spatial  infinity.   Let 

(6)  [y,u]   =  ^yaua  dt  dx  , 
and 

(7)  r[y]   =  <  ei[y'u]  >   =  f    ei[y'u]dP  . 


Then  F  is  a  functional  of  y  and  will  be  called  the  characteristic  functional 
of  P.   By  functional  differentiation  of  (7)  we  obtain 

;o\     "J /  .n  1/,     n     n/,     \  iLy,uJN 

— " =  <xu   (t1,x1)...u  (tn,xn)e      >  ; 

5y1(t1,Xl)-..5yn(tn,xn) 


hence 


(9)      (^(t^)...^^)  >  =  i"n- ^— 

Sy  ]  '   .   ...  .     • 


y=o 


The  left  side  of  (9)  is  just  the  n-point  space-time  correlation  function  of 
the  velocity  field. 

An  alternative  formulation  of  the  statistical  theory  of  turbulence  may 

be  based  on  a  sequence  of  probability  distributions  p  in  n-dimensional  spi 

n=l,2,...  .  These  distributions  are  defined  for  n  arbitrary  space-tj 

a      a 
(t, ,x  ),..., (t   ,x    )  and  n  arbitrary  components  u  , . . . ,u       of  the  vector 

a. 
field  u;  p(Sx, ...  j^Jd^, ...  ,dgn  being  the  joint  pro:  ■  i  '  J(t.,x.) 

are  in  the  intervals  (6.,|.+dg.);  j=l,...,n.   If  we  introduce  the  characteristic 

function 

r(o±,...,On)      =/exP   i^]  aj6j|p(fi1,...,fin)dg1,...,dSn 


(10 


of  this  probability  distribution,  then 


(11)         p(61,...,gn)    =    (2Jt)"n  Jexp  J-i  )t  o.z\  r(o1,...,on)&o1, 

et  us   evaluate  the  characteristic   functional  T  for  the  special 
argument  function  y  whose  components  are  given  by 


(12)  yP(t,x)      =     y^aJ5aJ3S(t  -tJ)5(x-xJ) 


Here  5   is  the  Kronecker  delta,  and  S(t)  and  5(x)  are  the  one -dimensional 

and  three-dimensional  Dirac  delta  functions.   For  y  given  by  (12), 

n 
(13)  [y,u]  =  Y]    "j^W' 


(lh)  T[y]   =  (expji^a.u  J(t.,x.)p-  T^, . .  .  ^   , 


and 


(15)     p(6-L,...,6n)  =  (2n 


/exPJ-i^]a.Jr[y]dai,...,dan. 


Thus  we  see  that  the  probability  distribution  functions  p  and  their  characteristic? 
functions  r   can  be  obtained  from  the  characteristic  functional  r[y]  by  means 
of  the  special  choice  (12)  of  y. 

3-   Derivation  of  the  functional  differential  equation. 

In  this  section,  we  obtain  equations  for  T   from  the  basic  flow  equations 
(2)  and  (3).   To  do  this  we  first  i   e  that 


/,/-\    of      /.  a,,    ,  i[y,u]  \        5T"  /  a,    >  0,,   ,  i[y,u]  \ 

(16)  =  (iu  (t,x)e  J'       )  ;    — = =  -  (u  (t,x)u  (t,x)e  lJ '        ), 

Sya(t,x)  8ya(t,x)5y^(t,x) 
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Hence,  from  (2), 


/,~\   5    5r      /  .  a  i[y,u]  \    /    .  r  Pa,     a     .a,   i[y,u]  \ 

(17)  ^7^77 —  =  (lv        >  =  <  l[-pa-uu  p  +  vuP  P+f  ]  e        > 

oy    (t,x)  x  x  xx 


a  =  1,2,3  • 
Now, 


3xV      Sya(t,x)  xV 


and 


o 

,.nv       5               5  r  /    r  a     0  ,    a  3   ,     i[y,u]  \ 

(19)     —a  — o =  -  (    [u  Qu    +u  u  J   e   L<y'       ) 

dxP   5ya(t,x)SyP(t,x)  xP                 xP 


But  the  quantity  u  „  in  (19)  vanishes  by  virtue  of  (3),  hence  the  last  three 


equations  yield 


(20)  I-      5r        =    i  -i ^ +  v  -j£ &L_  +  if°r  -  M.  . 

KeLK))   ot  _  a,.   *     X   .  P  _  a,.   w,  p,.   v  +   JJ  ,  a,t   \         *  a  ' 

By  (t,x)       dx   5y  (t,x)5y  (t,x)     dx  dx  &y  (t,x)  dx 


a  =  1,2,3 
n  =  <  ip(t,x)ei[y'u]  >. 


Various  devices  for  eliminating  the  pressure  term  -  — —  from  (20) 

dxa 

are  possible.   The  method  we  shall  use  here  is  to  introduce  the  "testing  field" 


r\(t,x)      =  (r\    ,T)  ,rf>) 


which  vanishes  sufficiently  rapidly  at  spatial  infinity  and  satisfies  the 
condition 


(22)  TfL   =  0  . 


(23 


)    f  f    iLta  -  -|,>=o 


and  from  (20)  we  obtain 


(2U)   n(t,x) 


a       &r        .    a  s2^ a2 


2 
3t   &ya(t,x)      X   BxP      6ya(t,x)5yP(t,x)       V  cbcPcbcP   &ya(t,x ) 

-  if^)  dt   dx     =     0   . 


This  is  a  functional  differential  equation  for  T[y].   It  must  be  satisfied  by 
r  for  all  testing  fields  r\   which  satisfy  (22).   An  alternative  to  the 
introduction  of  the  testing  field  r\   is  the  elimination  of  the  pressure  from 
the  Navier -Stokes  equation  at  the  outset.   By  using  the  boundary  conditions 
and  the  incompressibility  condition  .he  pressure  field  may  be  expressed  as 
the  sum  of  a  second-degree  polynomial  functional  of  the  velocity  field  and 
a  linear  functional  of  the  force  field.   (Cf.  [k] ,    for  example.) 

In  addition  to  (2k)   which  now  replaces  the  Navier -Stokes  equations  (2), 
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we  have  the  following  equation  which  replaces  the  incompressibility  conditior. 
(3): 


(25) 


cbc   &y  (t,x) 

This  is  an  immediate  consequence  of  (3)  and  (l6). 
Three  further  conditions  on  T   are 

(26)  r[o]    =    i,       r*[y]    =    r[-y]  ;  |r[y]|  <   1  . 

Here   denotes  the  complex  conjugate.   The  first  two  of  these  conditions 
follow  immediately  from  (7)  and  the  normalization  of  P.   The  third  condition 
is  a  consequence  of  the  normalization  and  non -negativity  of  P.   This  condition 
is  the  simplest  of  an  infinite  sequence  of  realizability  inequalities  which 
T  must  satisfy.   The  entire  set  of  inequalities  may  be  summarized  by  the 
statement  that  T[y]  must  have  an  everywhere  non -negative  functional  Fourier 
transform.   The  realizability  inequalities  imply  corresponding  inequalities 
for  the  velocity  moments. 

k.    Transformation  to  k -space. 

A  formal  simplification  of  the  functional  differential  equation  (2U) 
can  be  obtained  by  Fourier  transformation.   Let 

(27)  za(t,k)   =  jrya(t,x)e-ik'Xdx;   ya(t,x)   =  (2n)~3  J    ^(t^e^'^   ; 


(28)  Xa(t,k)   =  J   Tia(t,x)e"ik'X  dx  ;  ^(t,x)      =   (2tt)-3  f  >a(t,k)eik'X  dk 

(29)  ga(t,k)   =   (2Tr)-3^fa(t,x)elk'X  dx  ;  fa(t,x)  =  JV(t,k)e"ik-X  dk  . 


that 


(30)  z*(t,k)  =   z(t,-k)  ;   g*(t,k)  =  g(t,-k). 


An  immediate  consequence  of  (22)  is  the  condition 


(3D 


We  now  set 


(32)  rx[z]  =  r[y]  =  r[(2Tt)"3  ^    z(t,k)eik'Xdk]  . 

In  what  follows,  we  shall  omit  the  subscript  1,  since  it  is  always  clear  from 


(33) 


aPr 


5y1(Vx1)...5yn(tn,xn: 


■ — - —  exp  <  -i  (x,  •  k,  +  .  .  .  +  x  •  k  )  >  dk, 

ax  an         P    ll  1        n  n '(       1 


&z   (tn,kn  )...dz  "(t  ,k 
ll         n  n 


In  particular 


(3*>    -IT—    =    f 

Sya(t,x)  J       >z"(t,k) 


11 


e  dk 


5^  p  s2r  -ix-(k1+k. 


/ 


&ya(t,x)&yP(t,x)  J      ^(t^bt^it,^] 


e 


Now  each  term  of   (2k)  may  be  transformed  to  the  k-representation.      Thus 


(35)  J  t,    (t,x)  ^— — -dtdx  =   (2k)  Jj   a   (t,k»)e  gr  e  dk'dkdtdx 

^  5y    (t,x)  oz    (t,k) 


Txa(t,k)  £.-£ — at  dk 


5z"-(t,k) 


(36)    -v   /Y(t,x)   -|~3   — ^ dtdx 

J  BxV  ^ya(t,x) 


l(2rt)-3rxa(t^')eik'-XAP-/—    e^dk'dkdtdx  =vfr(t,k)k2-^- 
J  5za(t,k)  J  5za(t, 


— dtdk 
k) 


(37)      -i    /V(t,x)-\  5F       at  dx 

J  dxP  Sya(t,x)&yP(t,x) 


-ix-(k   +  k    ) 

dk, dk^dkdtdx 


(2rt)3^  X        2     6za(t,kl)&zP(t,k2) 


/>(t,k    +k    )(kP+k£)   — 2X_ dk  dkdt 

J  12        12      Sza-(tAi)5zP(t,k2)        1     2 


(38)  -i    T^f'Vdt  dx  ="I~  rxa(t,k')eik''Xga(t,k)e"ik'Xdk'dkdxdt 

=  -±t  rxa(t,k)ga(t,k)dtdk  =  -±r  rxa(t,k)[&aP-k"2kakp]g3(t,k)(itdk. 

The  last  equation  is  a  consequence  of  (31  )•   If  we  insert  (35  -  38)  in  (2*0  we 
obtain 

(39)  JVC*,*)  j  [|  ♦  ^"^  -  i[^-k-W]/(t,,)rJ  «  dk 


&2r 


-/At,k1+k2)(kJ+kp— -  -dkldk2d 

°  6z  (t,k1)&z  (t,k2) 
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This  is  the  functional  differential  equation  for  T[z].   It  must  be  satisfied 
by  T[z]  for  all  testing  fields  X  which  satisfy  (31  )•   In  the  k-space  repre- 
sentation, the  incompressibility  condition  (25)  becomes 

(ho)  f^F—    =    °> 

6z  (t,k) 
and  from  (26)  we  have  the  three  conditions 

(hi)  r[o]  =  i  ,    r*[z]  =r[-z]  ;   |r(2)|  <  1  . 

Our  explicit  introduction  of  the  projection  operator  6   -  k~  k  k  into  (39) 

emphasizes  the  fact  that  the  velocity  field  can  respond  only  to  the  transverse 
(i.e.,  solenoidal)  part  of  the  driving  forces. 

A  formal  representation  for  T[z]  is  the  "Taylor  expansion", 

(1+2)  T[z]  =  P  +  >  \  [t  (t..  ,k_  , . . .  ,t  ,k  ) 

v  o   / |  nJ  J      a  ,...,cl   v  1'  1*    *   n'   n; 

n=l 

n    a . 

v   I  |   z  J(t.,k.)dtndkn.  ..dt  dk 

X  I       v  j*  jy  1  1     n  n  ' 

j=l 

where 

c-np 

(1+3)  P        (tn,kn,...,t  ,k  N 


5Z  1(t1,k1 ),..., 6z  n(tn,kn: 


Now  from  (1+0)  we  see  that 


lU 


(kk)  kjJ  P^,...,^  =  °>  J=l,...,n; 


and  from  this  it  follows  that  there  exist  functions  Q       „   such  that 

1' '  '  '  '  n 


_nr       p  a   p  -, 
C*5)  P        (t^k,,...^  ,k  )  =TT  5„  R"k  k  k  K      r  (t,  ,k.,...,t  ,k  ) 


From  (Ul),  we  have  P  =  1  so  that  T  may  he  represented  in  the  form 

v*)  r  ^  =  i  +]T  sr/^  [\P;kv"2kXv]\,..vPn]T  ^^A' •  -dVkn  • 

n=l 
The  incompressihility  condition  (UO)  now  is  automatically  satisfied. 

From  (9)  and  (33)  we  see  that  the  space-time  correlation  function 

(hi)  U ■  (t1,x1,...,Vxn)  =(uai(t1,x1)...u\,Xn)  > 

1*    '  n 

is  given  by 

1'   '  n  J   1*  '    n 

For  many  purposes  it  is  useful  to  transform  (39)  so  that  a  characteristic 
Reynolds  number  appears  rather  than  tne  viscosity  coefficient  v.      This  can  be 
done  by  introducing  some  appropriate  characteristic  correlation  length  r  and 

characteristic  velocity  v  .   We  then  measure  lengths  and  velocities  in  units 

r^ 

This  last  unit  is  the  characteristic  time 


15 


If  we  denote  quantities  measured  in  the  original  units  by  primed  symbols  and 
corresponding  quantities  measured  in  the  new  units  by  un-primed  symbols,  then 

2 
r  v  v 

(k9)     t1  =  —  t  ,   x'  =  r  x  ,   u'  =  v  u  ,   k'  =  —  k  ,   f '  =  ~  f  , 


a-'   ,      v  a,       a'      v    a         a '   ,,.,      a3,  _. 

y  dx'  =  r-  y  dx  ,    z   =  p-  z   ,   and  g  dk'dt'  =  v  g  dk  dt 

v  r  v  r 

00  00 


From  (1+2)  and  (lj-3)  we  see  that 

&nr 


5z   , . . . ,&z    j=l 


T  z  ,:i(t.,k.)  dt.dk,  ...dt  dk 
II     K   3'   J         11     n  n 


is  dimensionless,  and  a  simple  calculation  shows  that  in  terms  of  the  new 
units  (39)  becomes 

(5o)  _/>(*.*)  {[I-2]^  -i%  -^W]/(t^)rj  at  a, 


dkxdk2dt  =  0. 

k2) 


Here  the  Reynolds  number  R  is  given  by 


v  r 
(51)         K  -  ■*? 
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5.  Solution  of  the  functional  differential  equation  for  zero  Reynolds  number 
If  we  take  the  limit  R  -> 0,  (50)  reduces  to 


(52)   /  At,k)<ls         --i|5nR-k-2kakV(t,k)rSdtdk 


/^^M^-^-^4 


In  this  case,  we  can  obtain  a  solution  of  the  functional  differential  equation 
and  recover  the  results  of  the  theory  of  "weak  turbulence",  derived  by 
Reissner  and  by  Batchelor  and  Townsend.    We  begin  with  the  following 

Lemma.   If  /  Xa(t,k)qa(t,k)dt  dk  =  0  for  all  X   which  satisfy  (31),  then 
there  exists  a  scalar  function  p(t,k)  such  that 

(53)  qa(t,k)   =  kap(t,k)  ;  a  =  1,2,3  . 


Proof: 


For  arbitrary   -r(t,k)  =  (t^t2,!-3),   set  Aa  =  ra  -  k"  k  t  ka 


^Ta(qa-  k-2kPqPka)dt  dk   =  JV"  k"2kPTPka)qadt  dk  =  J  x\a&t    dk  =  0  . 
Hence,  since  t  is  arbitrary,    q  =k(k  kq);a=  1,2,3  . 


If  we  apply  the  lemma  to  (52),  we  obtain 


^Mifc^v^-w-**' 


a  =  1,2,3 
'(t,k) 
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However,  from  (kO)   and  (5*0  we  see  that  pk  k  =  pk  =0.   Hence  p  =0  and  we 
obtain 

(55)  M  sfe "  l^v  k"w] gP(t,k)r  =  ° ;  a=1'2'3 ' 

From  (55)  we  see  that  satisfies  a  first  order  linear  inhomogeneous 

5za(t,k) 

ordinary  differential  equation  in  t  which  is  easily  solved  to  yield 


2 

(56)    — =  iha(t,k)r  +  e"k  l   Ba[z,k]  . 

52a(t,k) 


2    t  2 

(57)    ha(t,k)  =  (oap-  k-W)  e"k  t  f    gP(T,k)  ek  T  dr. 


and  from  (lj-0)  it  follows  that 

(58)  kaBa  =  0  . 
Now  the  substitution 

(59)  r[z]   =  Vtz]  exp  i±    /  ha(t,k)za(t,k)dt  dk> 

is  suggested  by  the  form  of  (56).   Indeed  (59)  is  the  analogue  of  the  usual 
exponential  form  of  solution  of  a  linear  ordinary  differential  equation,  and 
the  functional  V[z]  plays  the  role  of  the  coefficient  function  in  the  method 
of  "variation  of  parameters".   By  functional  differentiation  (59)  yields 
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(6o)        5F     =  iha(t,k)r  +  exp  i  /  hazadt  dk   - 

za(t,k)  J  Bz"(t,k 


and  by  comparing  (56)  and  (60)  we  see  that 

2 

(61)     — =  e     C  (z,k)  , 

&za(t,k) 

where 


=  exp  (-1    /  IvVdt  dk>  -a 


From  (62)  and  (58)  it  follows  that 

(63)  kaCa  =  0  . 

The  functional  C  is  arbitrary  so  far,  except  for  the  condition  (63)  and  the 
condition 


-k^   5^(z,kl)       -k|t2   &CP(z,k2) 
(64)   e       — p =  e       

&^(t2,k2)  5za(t1Jk1: 


which  follows  from  the  fact  that 


(65) 


6za(t1,k1)&zP(t2,k2)     SzP(t2,k2)6za(t1,ki: 
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Now  equation  (6l)  may  be  solved  for  V[z]  in  several  ways.   One  way 

a, 
is  to  expand  V  and  C  in  Taylor  series  and  compare  the  coefficients.   It 

is  easily  seen  that  the  solution  may  be  expressed  in  the  form 


(66)      V[z]   =  A[w];  wP(k)  =   f  e"kT(5   -k'2AP)za(T,k)dT 


(67)    5AO      =    f     e-«21     (&      K"2KaKP)  S(T-t)  6(K-k)dT 

5za(t,k)     J  aP 


2 
e'k  t(&ap-  k_2AP)  6(*-k)  , 


and  for  any  choice  of  the  functional  A[w], 


(68)   -^ =  [   _|A_  -^(JLL  dK  =  e"k2t(5  R-  k"2kV)  ~^- 

8za(t,k)     J      5wP(k)  6za(t,k)  aP  6w(i(k; 


Thus  (6l),  (63)  and  (61+)  are  satisfied  by  (66). 

We  now  insert  (66)  in  (59)  and  obtain  the  solution 

T[z]   =  A[w]  exp  \l    rha(t,k)za(t,k)dt  dk  i  ; 
(69)  [  J  J 

2   t         2 

,  cc, ,  ,  >    /c    ,  -2,  a,  0v  -k  t  f  P,  .  \  k  t, 
h  (t,k)  =  (&ap-k  k  kp)  e     /  g  (T,k)  e   dT  ; 


-  20  - 

2 
a./,  \      /<=    , -2n  a  J,   f  P,   ,  %   -k  t  , 
w  (k)  =   (S  B  -k  k  k  )   /  z  (T,k)  e     dT  . 

o 
The  functional  A  is,  so  far,  arbitrary  except  for  the  conditions 

(70)      A[0]   =  1,     A*[w]   =  A[-w],      |A(w)|  _<  1 

which  follow  from  (l+l)  and  (30).   It  is  easy  to  see  that  the  r[z]  given  by 
(69)  satisfies  (52)  and  (1+0). 

It  is  not  difficult  to  show  that  the  functional  A  is  simply  the 
characteristic  functional  (in  the  k-representation) of  the  distribution  of 
initial  velocity  fields  (cf.  Section  7).   This  distribution  must  be 
non-negative,  a  fact  which  implies  restrictions  on  the  possible  forms  of  A. 

For 
reduces  to 


(7D 


r[z]  -  a[w] 


functional  differentiation  we  obtain 


(72) 


6^ 


al 


5z"(tL^L)...SZ"(tn,kn) 

and  from  (1+3)  and  (1+8)  we  obtain 


ne  vv(wk  rVr 


S"A 


Sw     (k1)...Sw  "(k    ) 


73       U  (tn,x.,...,t  ,x 

v  '  ~"        an  ,...,av    1'    1'      '  n'    n 
1'       '    n 


n 


'•('.r^"'"1 


aB„...,p    <*i>-'V  exp^-K^-^  +...+  .Vkn)     dkr..dkn  . 


Here 


ilk)     aR      R  (k,,...,k  )  - 


5w  ±(k1)...5w  n(kn) 


Equation  (73)  provides  an  explicit  formula  for  the  time-dependence  of  the 
space-time  correlation  functions.   If,  for  example,  U       is  given  for 

t  =t  ,  ...,t  =  t  ,  then  aR      R  can  be  determined  by  inverting  the  resulti 

Fourier  transform. 

6.  Representation  of  averaged  Green's  function 

In  this  section  we  introduce  a  sequence  of  Green's  functions  which 
describe  the  averaged  response  of  the  velocity  fields  to  infinitesimal  pertur- 
bations in  the  external  force  field  f  .   These  functions  play  an  essential 
role  in  some  recent  studies  of  turbulence.   '   They  are  defined  as  follows: 


_  &u   (tnxn  ) 

GaP(Vx Jt2,x2)   ,     <  X1  > 

[75)  1122  s?»(t2x2) 


.  52ua(tn ,x.  ) 

GaPT(t     x,|t     x    -t     x    )      s    <    — ±-1 >      , 

112233  &fP(t   ,x    )SfT(t     x    ) 


Here  a, 3,...  =  1,2,3.  If>  as  we  have  assumed  in  section  2,  the  solutions 
u  (t,x)  are  determined  by  initial  conditions,  then  the  G  functions  vanish 
unless  t  >  t  ,  t  >  t  ,...  ;  that  is,  they  are  retarded  Green's  functions. 
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When  evaluated  at  f  =  0  the  Green' s  functions  give  the  averaged  response  of 
freely  decaying  turbulence  to  infinitesimal  external  perturbations. 

The  Green's  functions  may  be  represented  as  functional  derivatives, 
with  respect  to  y  and  f,    of  the  solution  T   of  {2k).      Noting  (l6),  we  have 


(75a) 


cf3(vXl|t, 


&fP(t2,x2) 


£T 


&JT  (t1,X1) 


y=o 


Gaf3r(trx1|t2,x2;t3,x3: 


.-1 


6fp(t0,x0)6f'(t,xJ 


3'  3' 


or 


By  (t1,x1: 


y=o 


In  the  "weak  turbulence"  limit,  R  -»  0,  we  may  use  the  results  of  the 
preceding  section  to  calculate  the  Green's  functions.   Only  the  lowest  member 
G   of  the  sequence  now  is  non- vanishing.   The  higher  Green's  functions  describe 
intrinsically  non-linear  effects  which  are  ignored  in  this  approximation. 
In  order  to  calculate  G   we  combine  (69)  and  (29)  to  obtain 

(76)   ha(t1,k1)  =  (2*r3(oap-  k"2k^)  y^cLt  ^dx  exp  Jk^(t-t1)+ik1-xi  fP(t,x)  , 


and  functional  differentiation  of  this  equation  yields 
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61^  ,kn  ) 


(77)       p     1?   1     =   (2n)-3(^-k'12k<^)J dt  J dx  exp  <kj(t  -t^  +  11^x1  Bfc-t2)6(x  -x2) 


6fp(t2,x2) 


(2*r3(5  -k^2k^)  exp  ^(tg-^)  +l*1-x2>  if  t2  <  tx| 


Now  from  (69)  we  see  that 


(78) 


Sz^W 


-k2tn 
ih^t^k,  )  +  e  L  L 


l'T. 


P\  &A 


tp  kl  klkl  I 


&wp(kx) 


hence 


(79) 


&fH(t2,x2)  l&zu'(t1,k1) 


or 


&hu'(t1,k1) 


It  follows  that 


(80)  G 


1  ll  2  2      ^(t2,x2)J         K<Wz-oJ 

r  -ixi-ki  5"a(ti'ki'  .. 


-  2k  - 


3  /  f5ap-kI2klk?)   eXP  \ -k?(tl't2 }  +  ikl' (X2  "  Xl } \   dkl    if  *2  -  *! 


7.  Treatment  of  random  driving  forces 

We  wish  now  to  describe  briefly  the  generalization  of  the  functional 
formalism  to  the  case  where  the  driving-force  field,  as  well  as  the  velocity 
field  is  statistically  distributed.   We  shall  discuss  here  only  the  case  where 
the  force- field  and  the  initial  velocity  field  are  statistically  independent. 
Instead  of  T,  we  consider  the  joint  characteristic  functional 


(81)   H[y,q] 


i[y,u]  +  i[q,f]    > 


where  y(t,x)  and  q(t,x)  are  arbitrary  vector  fields  and  the  averaging  <   )  now 
is  over  the  joint  probability  distribution  of  u  and  f.   The  characteristic 
functional^  of  the  initial  velocity  distribution  and  of  the  force  field  distribution 
are  defined,  respectively,  by 


(82)      A[v]   =  <  exp[i  /  v(x)  u(0,x)  dx]  ) 


and 


(83)      *[q]   =  <   ei[^f]   > 


Here  v  is  an  arbitrary  vector  function  of  x.   We  regard  A  and  \\r   as  prescribed. 
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It  follows  immediately  from  the  statistical  independence  of  force  field  and 
initial  velocity  field  that  ft  must  satisfy  the  "boundary"  condition: 


(81+)    ft[y,q]   =  A[v]   i|r[q] 


If     y(t,x)  =  5(t  -0+)v(x)    . 


In  particular, 

(85)        ft[0,q]      =     *[q]      . 

All  moments  of  the  joint  distribution  of  velocity  and  forcing  fields 
may  be  expressed  as  functional  derivatives  of  ft.   For  example,  we  have 


(86)   <  fa(t,x)  ) 


.-1 


&q  (t,x) 


y,q=o 


<  ua(t,x)  fP(t',x' )  >  =  i" 


Oy"(t,x)SqM(t  ,x 


y,q=o 


If  now  we  retrace  the  derivation  of   (2k)  and   (25),   we  find 


(87)    jY(t,x) 


axpaxpy    5ya(t,x) 


5ft 


m 


axP   5ya(t,x)&yP(t,x)        &qa(t,x) 


dt  dx     =     0 
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and 


(88) 


cfac   Sy  (t,x) 


It  is  desired  to  find  the  solution  of  (87)  and  (88)  which  satisfies  (8k). 

In  order  to  illustrate  the  formalism,  we  shall  give  the  solution  in 
the  limiting  case  of  vanishing  Reynolds  number.   (We  choose  to  retain  the 
x-space  representation  here,  rather  than  make  the  transformation  to  the 
equivalent  k-space  formalism. )  In  the  zero  Reynolds  number  limit,  (87) 
becomes 


JVtt.x)  I 


&-th)  PZrr^ 


where  the  units  introduced  in  (k9)   are  used.   The  solution  of  (89)  which 
satisfies  (8k)   and  (88)  is 

(90)    «[y,q]   =  A[r]  ^[q+r]  , 


where 


(92 


)   ra(t,x)  =   /  yP(t',x' )GPa(t',x' |t,x)dt'  dx' 


-  27 


The  validity  of  (90)  may  be  verified  by  substitution  and  use  of  the  definitions 
(82)  and  (83).   We  shall  only  outline  the  demonstration  here.   Equation  (8k), 
first  of  all,  follows  directly  from  the  expression  (80)  for  the  Green's  function. 
The  various  functional  derivatives  which  occur  in  (88)  and  (89)  may  be  expressed 
in  terms  of  f  and  G  by  use  of  the  chain  rule  for  functional  differentiation. 
Equation  (88)  then  follows  immediately  from  (80 ).   Equation  (89)  follows  from 
(80)  and  (22),  after  an  appropriate  partial  integration  is  performed. 

8.  Conditions  for  homogeneity 

We  wish  finally  to  see  how  the  special  conditions  of  homogeneous  turbulence 
appear  in  the  present  formalism.   Homogeneity  is  defined  in  terms  of  the 

probability  distribution  P  introduced  in  section  2.   P  is  said  to  be  homogeneous 

12  3 
if  it  is  invariant  under  all  translations  L  ;  a  =  (a  ,a  ,a  );   where 

(93)  L  u(t,x)   -  u(t,x  +  a)  . 

From  (6)  and  (7)  we  see  that  P  is  homogeneous  if  and  only  if 

(9k)  T[Lay]   =  r[y]  ,        for  all  a; 

and,  of  course,  this  condition  is  equivalent  to  the  condition 


(95)  r[L  z]   =  r[z]  ,       for  all  a; 


where 


(96) 


f  /      \  -ik«x,      r  ,.   v  -ik- (x-  a)   ia-k  ,.  ,  ^ 
L  z(t,k)  =  /  y(t,x+a)  e     dx  =   /  y(t,x)e         =  e    z(t,k). 
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From   (U6)  we   see  that 


n=l  v=l 


Oo,     P 

a  P       v     v      v 
v  v 


n       ia-k .  a. 

,    e  Jz  J  dt_.dk...  dt 

P-,,...,P     M  linn 

J=l 


v-aH 


Thus  (95)  implies  that  for  all  a  and  n=1,2,... 


^8)  \,...,^\>\>"> W  =  6ZP  ^-(V-^j  \,...,pfl^— Vkn: 


It  follows  that  another  equivalent  condition  for  homogeneity  is  that  there 
exist  functions  HR     R  (t  ,k ,,...., t  ,k  )  such  that 


(99)  Q, 


^...^W'^nV  =   5(V---+kn}  V..,P^kl-tn'kn: 
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